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ABSTRACT 

We have obtained long-slit spectroscopy for a sample of 9 SO galaxies in the 
Fornax Cluster using the F0RS2 spectrograph at the 8.2m ESO VLT. From these data, 
we have extracted the kinematic parameters, comprising the mean velocity, velocity 
dispersion and higher-moment and /14 coefficients, as a function of position along 
the major axes of these galaxies. Comparison with published kinematics indicates that 
earlier data are often limited by their lower signal-to-noise ratio and relatively poor 
spectral resolution. The greater depth and higher dispersion of the new data mean that 
we reach well beyond the bulges of these systems, probing their disk kinematics in some 
detail for the first time. Qualitative inspection of the results for individual galaxies 
shows that they arc not entirely simple systems, perhaps indicating a turbulent past. 
Nonetheless, we are able to derive reliable circular velocities for most of these systems, 
which points the way toward a study of their TuUy-Fisher relation. This study, along 
with an analysis of the stellar populations of these systems out to large galactocentric 
distances, will form the bases of future papers exploiting these new high-quality data, 
hopefully shedding new light on the evolutionary history of these systems. 
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1 INTRODUCTION 

One of the key areas of research in extragalactic astron- 
omy is the study of the formation and evolution of galaxies 
in different environments, from the low-density field to rich 
clusters. Our understanding of the individual physical pro- 
cesses involved in the evolutionary history of galaxies, and 
their relative importance in each environment, is still rather 
poor. In this context, lenticular (SO) galaxies have a very 
important role. Located between elliptical and spiral galax- 
ies in the Hubble Diagram, SO galaxies have been a focus of 
debate for many years. One fundamental - and sometimes 
controversial - issue is whether the formation of these galax- 
ies is more closely linked to that of ellipticals or to that of 
spirals. 

The evolution with redshift of the morphology-density 
relation in clusters of galaxies (Dressier 1980; Dressier et 
al. 1997) implies that the relative fraction of SO galaxies 
increased from z ~ 0.5 to the present, while the relative 
fraction of spirals decreased in a similar proportion. An 
evolutionary connection between SOs and spirals has thus 
been proposed, at least in cluster environments. More re- 
cent works seem to support such ideas (e.g., Fasano et al. 
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2000; Desai et al. 2006), although alternative views do exist 
(Andreon 1998). However, our understanding of the mecha- 
nisms involved in the possible morphological transformation 
of spirals into SOs remains poor. In addition, close similari- 
ties exist between ellipticals and SOs in their colours, stellar 
populations, gas content, location on the fundamental plane, 
etc, so the debate as to whether SOs are more closely related 
to spirals or ellipticals remains open. 

The presence of SOs in both cluster and field environ- 
ments raises the real possibility that multiple evolutionary 
paths exist for the formation of these systems. Indeed, a vari- 
ety of mechanisms have been proposed that can, in principle, 
alter a galaxy's morphology in this direction. These mech- 
anisms include ram pressure gas stripping (Gunn & Gott 
1972), harassment (Moore et al. 1998), starvation (Larson, 
Tinsley & CaldweU 1980; Bekki et al. 2002), unequal mass 
galaxy mergers (Mihos & Hernquist 1994; Bekki 1998) and 
galaxy-cluster interactions (Goto et al. 2003). However, it is 
still unclear whether any of these processes would work in 
practice, and, if more than one, their relative importance in 
different environments. 

Van den Bergh (1990), analysing the Revised Shapley- 
Ames Catalog of Bright Galaxies (Sandage & Tammann 
1987), found that the frequency distribution of the lumi- 
nosity of SOs is not intermediate between E and Sa galax- 
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ies. This discontinuity could imply the existence of sub- 
populations amongst the SOs: bright SOs, the "real" interme- 
diate class between E and Sa galaxies; and faint SOs, many 
of which could be miss-classified as faint Es if viewed close 
to facc-on. The works of Nioto ct al. (1991) and Jorgcnscn & 
Franx (1994) support this hypothesis, pointing out the sim- 
ilarity between disky [and thus faint (Kissler-Patig 1997)] 
E and SO galaxies, based on their isophotal central shapes. 
Also, Graham et al. (1998), in a study of the extended stel- 
lar kinematics of elliptical galaxies in the Fornax Cluster, 
found that five of the galaxies are in fact rotationally sup- 
ported systems, suggesting that they could be misclassified 
SO galaxies. 

Studies of the stellar populations in cluster galaxies 
(Kuntschner 2000, Small et al. 2001) also support the idea of 
a dichotomy between low- and high-luminosity SOs: bright 
SOs are old and coeval with E galaxies, while faint members 
present younger central ages, indicating more recent star 
formation episodes. Furthermore, Poggianti ct al. (2001) ex- 
amined the star-formation history of early-type galaxies in 
the Coma Cluster, and found that ~ 40% of the SO pop- 
ulation seemed to have experienced a star-formation event 
during the last few billion years, a phenomenon which is 
absent in their sample of elliptical galaxies. Thus, it has 
been proposed that faint SOs could be the descendants of 
the post-starburst galaxies found in high redshift clusters. 
The work of Mehlert et al. (2000, 2003) in early type galax- 
ies in the Coma Cluster confirms the dichotomy found by 
Poggianti et al. between old and young lenticulars; however, 
the high alpha-element ratios found in the latter seems to 
argue against the occurrence of recent star formation: the 
authors suggest that the strong Balmer line indices mea- 
sured in apparently "young" SOs could actually be produced 
by unusually blue horizontal branches rather than by young 
stellar populations. Clearly, more work remains to be done 
in the study of stellar populations to interpret the physical 
significance of the two apparently-distinct types of SO. 

The formation of lenticular galaxies has also been the 
subject of numerical simulations. The work of Shioya et al. 
(2004) on "red H(5-strong" galaxies also suggests two differ- 
ent evolutionary paths for SOs, each one able to match dif- 
ferent spectral features in different galaxies: a "truncation" 
scenario, in which the star formation is stopped and followed 
by passive evolution, and a "starburst" one, in which a rel- 
atively recent and short star formation event precedes the 
cessation of star formation. On the other hand, Christlein & 
Zabludofi^ (2004) found that simulations based on a fading 
stellar population could not match the observed luminos- 
ity distribution of galaxies with the larger bulge-to-total- 
light ratios {B/T) typical of SOs. They therefore advocated 
a "bulge enhancement" model, where disk fading is accom- 
panied by an increase in the luminosity of the bulge; this 
model seems to match the observations over a wide range 
of B/T. Although this work still faces the problem of the 
difficulty of reliably identifying SO galaxies, it does raise the 
interesting possibility that, at least in some cases, the cessar- 
tion of star formation in the disk of the parent spiral galaxy 
could be preceded by a centrally-concentrated episode of 
star formation, thus enhancing the bulge luminosity. In this 
context, it is interesting that Moss & Whittle (2000) found 
spiral galaxies in clusters have more central star formation 
than their field counterparts. 



The extended stellar kinematics of SO galaxies have 
been studied mostly in conjunction with those of ellipti- 
cals, and often as one single class of objects, generically 
termed "early-type galaxies." For instance, D'Onofrio et al. 
(1995) studied a sample of 15 early-type galaxies in the For- 
nax Cluster, and did not find major differences between Es 
and SOs other than stronger rotational support and higher 
projected ellipticities for the latter. In a sample of 35 E 
and SO galaxies in the Coma cluster, Mehlert et al. (2000, 
2003) found that elliptical galaxies have, on average, slightly 
higher velocity dispersions than SOs, as is also apparent from 
the velocity dispersion profiles presented in D'Onofrio et al. 
(1995). Although these differences could be real, they may 
also at least in part arise from the selection effects that ren- 
der SOs more reliably identified when close to edge-on. 

With current techniques using integral-field units and 
high quality spectra, it is possible to examine stellar- 
kinematic substructure in search of further clues as to how 
these systems form. The work by Emsellem et al. (2004), 
for example, revealed that kinematically-decoupled compo- 
nents, bars and misalignments between photometric and 
kinematic axes seems to be present in both Es and SOs. 
There arc also instances of even more extreme kinematic 
substructure such as counter-rotating co-spatial stellar disks 
in SOs (Rubin, Graham & Kenney 1992), but these seem 
to be very rare (Merrifield, Fisher & Kuijken 1996). This 
rarity is something of a surprise, as counter-rotating gas is 
relatively common in SOs (Bertola, Buson & Zeilinger 1992), 
which led to the suggestion that it might be quite common 
for SOs to be enhanced by the kind of minor mergers likely re- 
sponsible for this phenomenon. These observations can only 
be reconciled if there is some mechanism for inhibiting star 
formation in such counter-rotating material, but the situa/- 
tion is clearly quite complex. 

One last approach to understanding the origins of SO 
galaxies lies in their scaling relations. In particular, there 
are a handful of studies of the TuUy-Fisher relation for SO 
galaxies. Using a sample of field SO galaxies, Neistein et al. 
(1999) found a much larger scatter in this relation between 
luminosity and rotation speed than is found in spirals. Hinz 
ct al. (2001, 2003) found a similar result for cluster SOs. The 
large scatter would tend to indicate a rather stochastic his- 
tory for these systems rather than any simple deterministic 
evolutionary path from spirals. By contrast, Mathieu, Mer- 
rifield & Kuijken (2002) selected a sample of SOs with rel- 
atively small bulges and found a TuUy Fisher relation with 
a small scatter but an offset to fainter luminosities relative 
to the spiral relation, suggesting a relatively simple fading 
of spirals into SOs for this sub-sample of systems. 

The above discussion illustrates the wide range of tech- 
niques that have been applied to trying to understand the 
nature of SOs, and the rather contradictory results that have 
emerged. There are certainly indications of a dichotomy be- 
tween faint and bright SOs, but it is still unclear which ob- 
servables best characterize this distinction, and how those 
observables might translate into differences in the evolution- 
ary history of the two types. Much of the difficulty arises 
from the heterogeneous nature of the data that have been 
used in these studies. In some cases, the data come from ob- 
jects in a range of ill-defined differing environments, while 
in others the issues are more to do with the varying quality 
of the observations. 
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To address these issues, we have obtained very high 
signal-to-noise ratio long-slit spectroscopy of 9 SO galaxies 
in the Fornax Cluster. The data are uniform in quality and 
probe in detail a single environment, removing these vari- 
ables from the problem. However, the sample does include 
both faint and bright examples of SOs, allowing us to investi- 
gate whether there are systematic differences between these 
systems in anything other than luminosity, and hence ad- 
dress the question of whether they formed in different ways. 
In this paper, we present the details of the sample, the obser- 
vations, the data reduction, the resulting kinematics, and a 
qualitative look at the results. In subsequent papers, we will 
delve deeper into the quantitative information that these 
data provide, looking at their TuUy-Fisher relation and the 
nature of their stellar populations out to large radii. 

The layout of this paper is as follows: Section |21 de- 
scribes the sample and observations; Section |3 explains the 
basic data reduction; Section 2] describes the extraction of 
the kinematics; Section |S] presents the resulting kinemat- 
ics and comparison to previous data; Section |S| extracts the 
rotation curves from these data; and Section Q contains a 
summary. 



2 THE SAMPLE AND OBSERVATIONS 

The sample was selected from galaxies in the Fornax Clus- 
ter classified as SOs by the NASA/IPAC Extragalactic 
Database. They were selected to span a wide range of lu- 
minosities ( — 22.3 < Mb < —17.3), and to be sufficiently 
inclined to measure rotations along their major axes. The 
basic properties of the resulting sample of SOs are presented 
in Table 1. 

The necessary observations of the major axes of these 
galaxies were made in service mode between 2002 October 
2 and 2003 February 24 at the 8.2m Antu/VLT using the 
F0RS2 instrument in long-slit spectroscopy mode; exposure 
times and dates are provided in Table 1. Spectrophotomet- 
ric standard stars were observed each night, and we also 
obtained spectra of stars with a range of spectral classes to 
act as velocity templates in the kinematic analysis; these 
objects are listed in Table 2. During the observations, the 
seeing varied from 0.75 to 1.48 arcsec FWHM, which is more 
than adequate for the study of these large objects. 

The detector in F0RS2 comprises two 2kx4k MIT 
CCDs, with a pixel size of 15 x 15/um^. The standard res- 
olution collimator and the unbinned readout mode were 
used, yielding a scale of 0.125"/pixel. The spectrograph 
slit was set to 0.5" wide and covered 6.8' in length. The 
GRIS1400V -I- 18 grism was used, providing a dispersion of 
0.318A/pixel and covering the 4560A A < 5860A wave- 
length range. This set-up provided a spectral resolution, as 
measured from the FWHM of the arc lines, of « 4 pixels 
(or 1.12 A), which translates into a velocity resolution of 
73.3 kms"^ FWHM (or 31.0 kms"^ in terms of the veloc- 
ity dispersion). The CCD was read out at 200kHz, which 
is twice the normal speed used for spectroscopy. The high 
readout speed was the only one available for unbinned read- 
out of the CCD. 



3 DATA REDUCTION 

The standard data reduction process was carried out using 
IRAF (Tody 1986, 1993). Bias subtraction, flat-fielding and 
cosmic ray removal were applied to all science and calibra- 
tion images. Bad pixel masks were created by dividing two 
fiatfield images of different exposure times, and the bad pixel 
values were interpolated. When more than one exposure was 
obtained for each galaxy (see Table 1), the individual spec- 
tral images were combined to maximise the signal-to-noise 
ratio (S/N). During the combination process the position of 
the galaxy spectrum and the locations and widths of sev- 
eral sky lines were checked; the match between the different 
exposures for each galaxy was found to be excellent, so no 
further alignment was necessary. 

After tracing the arc lines of the He-Ne lamp, a 2D 
correction was applied to all the spectra in order to remove 
the geometric distortions due to the instrument optics and 
CCD flatness issues. The wavelength calibration was applied 
in the same step. The residuals of the wavelength fits were 
typically 0.1-0.2A. 

Next, we subtracted the night sky spectrum using the 
IRAF task background. Potentially, residuals and uncertain- 
ties from sky subtraction could bias the kinematic and line 
strength measurements that we will make from these data, 
so we have looked quite closely at how this step is imple- 
mented. Given the large spatial coverage of the data, there 
are sizeable regions that are free of galaxy light from which 
we can define a sky spectrum, and we found that the result- 
ing sky-subtracted galaxy spectrum was almost completely 
insensitive to how the sky region was chosen. We have there- 
fore opted for the simplest choice of selecting sky regions, 
typically 100 pixels wide, on each side of the galaxy and 
as far from it as possible within the confines of the slit. 
The only exception is the case of NGC 1316, the largest and 
brightest galaxy, where only one sky region was used because 
the galaxy centre was placed close to one end of the slit in 
order to maximise the spatial coverage. Contamination by 
scattered light was not found to be a significant issue: in 
all cases, there were no nearby bright sources, and the rel- 
atively low surface brightnesses of the systems reduces the 
likelihood of misplaced light from the galaxies themselves. 

A periodic, square wave noise pattern or 'hum' was 
found to be present along the dispersion direction in all 
the spectra. The amplitude of this hum varied from 3 to 
6 counts, and its frequency, 0.077 cycles/pixel (15kHz), was 
reasonably constant, with small deviations due mainly to 
previous data processing. Because this pattern could par- 
tially mask the spectral features at low S/N, we subtracted 
it using a Fourier technique (Brown 2001). Ideally, one would 
apply this technique to the raw data frames, but its low level 
meant that it could only be reliable identified and removed 
after processing to subtract the sky. In practice, inspection 
of the resulting cleaned images showed that this processing 
had no adverse impact on our ability to remove this source 
of systematic noise. 

To correct by the effect of atmospheric extinction, we 
used the IRAF task setairmass to define the effective air- 
mass for each spectrum. Then, we applied the extinction cor- 
rection (task calibrate) using the appropriate atmospheric 
extinction table. The spectra were subsequently fiux cali- 
brated using the set of spectrophotometric standard stars 
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Table 1. Sample of SO galaxies in Fornax. 



Name RA DEC Diameter* Exp. Time Date 

['] [sec] 



NGC 1316 


03 


22 


41 


-37 


12 


30 


9.40 


11.0 


3 


X 


1200 


13 


Oct 


2002 


NGC 1380 


03 


36 


27 


-34 


58 


34 


10.9 


4.8 


2 


X 


1200 


24 


Feb 


2003 


NGC 1381 


03 


36 


31 


-35 


17 


43 


12.4 


2.7 


2 


X 


1600 


24 


Feb 


2003 


IC 1963 


03 


35 


30 


-34 


26 


51 


12.9 


2.6 


2 


X 


1600 


31 


Jan 


2003 


NGC 1375 


03 


35 


16 


-35 


15 


56 


13.2 


2.2 


2 


X 


1800 


28 


Dec 


2002 


NGC 1380A 


03 


36 


47 


-34 


44 


23 


13.3 


2.4 


2 


X 


1700 


28 


Dec 


2002 


ESO358-G006 


03 


27 


18 


-34 


31 


35 


13.9 


1.2 


2 


X 


2400 


14 


Oct 


2002 


ESO358-G059 


03 


45 


03 


-35 


58 


22 


14.0 


1.0 


1 


X 


2550 


8 


Feb : 


2003 


ESO359-G002 


03 


50 


36 


-35 


54 


34 


14.2 


1.3 


1 


X 


2250 


26 


Nov 


2002 



* From RC3, de Vaucouleurs et al. (1991). 



Table 2. Spectrophotometric (S) and Template (T) Stars. 



Name 


T/S 


Spectral Class 


BD-01-0306 


T 


Gl V 


HD1461 


T 


GO IV 


HD7565 


T 


K2 III-IV 


HD16784 


T 


GO V 


HD18234 


T 


KO 


HD19170 


T 


K2 


HD21197 


T 


K5 V 


HD23261 


T 


G5 


HD36395 


T 


Ml. 5 V 


HD61606 


T 


K2 V 


LTT1020 


S 




LTT1788 


S 




LTT3218 


S 




HZ4 


S 




GD71 


s 




Feige67 


s 





listed in Tabic 2, which allowed us to transform the ob- 
served counts into relative spectral fluxes as a function of 
wavelength. The sensitivity functions as derived from the 
different stars varied by typically ~ 1%, reaching ~ 3% to- 
wards the edges of the wavelength range, so we created a 
single function from all of them and can be reasonably con- 
fident that the resulting internal relative flux calibration is 
good to ~ 1%. 

This analysis was applied both to the galaxy spectra 
and to the spectra of the 10 template stars (listed in Ta- 
ble 2), which will be used to model the spectra of the galax- 
ies during the extraction of the kinematics. 

The calibrated two-dimensional galaxy spectra were 
binned along the spatial direction so as to generate one- 
dimensional spectra from different radii with comparable 
S/Ns. We started at the centre of the each galaxy using bins 
with a S/N per Angstrom of 53 (or 30 per pixel). When we 
reached the outer parts, where no more bins with that S/N 
ratio could be built, we reduced the S/N of each bin to 30 
per Angstrom, then to 15, and finally to 10. The data were 
also re-binned on to a logarithmic wavelength scale in the 
dispersion direction to enable the kinematic measurements. 



4 EXTRACTION OF THE KINEMATICS 

In the present study, the kinematic properties of the galax- 
ies were derived using the Penalised Pixel Fitting method 
(pPXF, Cappellari & Emsellem 2004), which is a parametric 
technique based on maximum penalised likelihood. Since the 
resulting kinematic parameters are so central to this work, 
we first describe the workings of this method as they apply 
to the current data set, and the tests that we undertook to 
check the reliability of the resulting parameters. 

The parametric method used by pPXF models the line- 
of-sight velocity distribution (LOSVD) as a Gaussian plus 
a series of Gauss-Hermite polynomials. The fits return the 
mean line-of-sight velocity, VLos, and the velocity disper- 
sion, a, from the Gaussian, together with the /13 and /14 
coefficients from the polynomials. These two coefficients are 
related, respectively, to the skewness and the kurtosis of 
the LOSVD, which are higher moments associated with the 
asymmetric and symmetric departures of the LOSVD from 
a Gaussian. 

This software works in pixel space, finding the com- 
bination of template spectra which, convolved with an ap- 
propriate LOSVD, best reproduces the galaxy spectrum in 
each bin. A subroutine of the program fits the continuum, 
using Legendre polynomials of order 4 (in this case), and 
divides the spectra by the fit to remove any possible low- 
frequency mismatches between the galaxy spectra and the 
model. The method carries out a pixel-by-pixel minimisation 
of the residuals (quantified by x^)> adding an adjustable 
penalty term to the x'^ to bias the resulting LOSVD to- 
wards a Gaussian shape when the higher moments /13 and 
/14 are unconstrained by the data. Thus, a deviation D of the 
LOSVD from a Gaussian shape will be accepted as an im- 
provement of the fit only if it is able to decrease the residuals 
by am amount 

Xl = x\l + X'-D'). (1) 

The A parameter defines the threshold between a deviation 
considered as an improvement of the pure Gaussian fit, and 
one considered as due to noise, and thus rejected. As Cap- 
pellari & Emsellem (2004) suggest, the value of A has to be 
determined for each particular dateiset because the ability 
of the software to find the correct LOSVD parameters (es- 
pecially the higher moments) will depend critically on the 
S/N of the data and other properties of the particular spec- 
tra such as the spectral range covered. To estimate the op- 
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Figure 1. Sample 1000 pPXF simulations carried out in order to 
find the optimum value of the A parameter for this dataset. In this 
particular case, the simulated spectra have S/N = 53 and A has 
been set to 0.3. The model LOSVD has Vi^ = 0, ftsin = -0.07, 
^14111 = 0.09 and the derived values of Vout, hzo\it and /i4out are 
shown for different values of (Tin. 



timum parameter for our study, we performed Monte-Carlo 
simulations in which we convolved the template spectra with 
a LOSVD whose shape was specified by model parameters 
Vin, /i3in and /i4in, and whose width was allowed to vary 
such that < CTin < 250 kms^^. These spectra, degraded to 
the S/N of the various one-dimensional galaxy spectra, were 
then passed through the pPXF software to recover estimates 
for the parameters, Vout, fout, /isout and ^4out. An example 
of one such run is presented in Figure^ After extensive test- 
ing using a variety of input parameters and S/N, we found 
that A = 0.3 provided an optimum value for the spectra 
in the current dataset, yielding values that minimized both 
errors and bias in the results. 

As is apparent from Figure even at this fairly modest 
value for the A parameter, a degree of bias appears in the 
answers at low dispersion, particularly in the derived values 
of cr and hi. This kind of bias, particularly prevalent in such 
even moments, is a well-known phenomenon (e.g. Cappellari 
& Emsellem 2004). Clearly, we must understand and possi- 
bly correct for this effect before interpreting the calculated 
parameters, so we have performed a range of tests using our 
adopted value of A but allowing other parameters to vary. 
The conclusion of this study is that the bias in a is small 
but well-defined at a maximum level of ±2-3kms~^; we cor- 
rected the results for this bias, but it makes no substantive 
difference. The case for hi is a little more disturbing, as this 
parameter, particularly in galaxies with low values of a, can 
suffer significant bias, amounting to offsets as high as 0.4 in 
the most extreme cases. The uncertain nature of this bias 
means that we have not corrected for it explicitly in our fi- 
nal derived value, but we would caution against interpreting 
the hi values as too precise, particularly from spectra where 
S/N is low and where a is small. Nonetheless, this study 
does demonstrate that the principal parameters of interest, 
Vlos and a, are robustly extracted from spectra of the qual- 
ity that we present here, irrespective of the detailed shape 
of the LOSVD as specified by /13 and hi. 



To estimate the errors on the derived parameters for 
each one-dimensional galaxy spectrum, we ran a further se- 
ries of Monte Carlo simulations. We modelled each such 
spectrum using the corresponding best combination of tem- 
plate stars, broadened by the best-fit LOSVD, and degraded 
to the S/N of the galaxy data. These model spectra were 
then passed through the pPXF software to measure the dis- 
tribution of the derived parameters, and hence a measure of 
the uncertainty in a single spectrum. The dispersion in the 
parameters derived in each Monte Carlo simulation is quoted 
as the error bar on the corresponding parameter derived for 
the real galaxy spectrum. 

Before we can present the final data set of kinematic 
parameters versus distance along the major axes for these SO 
galaxies, the only remaining issue is to define the zero-points 
of velocity and position for each galaxy. The velocity zero- 
point was simply set to the median value of VLos determined 
along the slit for each galaxy. For the spatial zero-point - 
the galaxy's centre - the simplest definition would just be 
the point along the slit that admitted the greatest amount 
of light. However, this definition could well pick up on an off- 
centre localized feature that would then distort the over-all 
kinematic structure of each galaxy that is the primary goal 
of this study. Instead, we adopted as the spatial origin the 
position along the slit that renders the plot of Vlos versus 
radius as close to anti-symmetric as possible. The only case 
where this method was not applied was NGC 1316; the large 
size of this system meant that we placed the slit primarily to 
cover only one side of the galaxy, so we could not subtract 
the median velocity along the slit nor match up the two 
sides of the rotation curve to find the kinematic centre. For 
this galaxy, we used the point along the slit at which we 
derived the largest value of cr to define the spatial location 
of the galaxy's centre, and then used the value of Vlos at 
this point to specify the galaxy's systemic velocity. 



5 RESULTS AND COMPARISONS 

The values for the kinematic parameters Vlos, o", hz and 
hi versus radius, as derived using the analysis set out in 
Section^ are presented in Appendix El 

A number of these galaxies have been studied previ- 
ously with lower quality spectra or less sophisticated spec- 
tral analyses, so as well as presenting our own results we 
will also compare the data to these published results. The 
most extensive study for comparison is that of D'Onofrio et 
al. (1995), who obtained major-axis kinematics (Vlos and 
cr) for six galaxies in common with the current sample. One 
complication in comparing their data with the new obser- 
vations is that they presented the kinematics folded about 
the centre of each galaxy, leading to an ambiguity in the 
actual spatial location of each measurement. We therefore 
cannot match up individual features that occur on one side 
of a galaxy, but can use these data to examine the overall 
profile. The other limitation of their data was its relatively 
low spectral resolution (~ 60kms~^ in velocity dispersion), 
which, as we will see, introduces a significant bias into the 
kinematics. Further long-slit data for single galaxies in the 
current sample can be found in Graham et al. (1998) and 
Longhetti et al. (1998), so we also include these data in our 
comparison. There are also a number of measurements of the 
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central velocity dispersion that provide a single comparison 
point. Kuntschner (2000) obtained central dispersions for all 
the galaxies in this sample, although he warns that the rel- 
atively low spectral resolution of his data, corresponding to 
a velocity dispersion of ~ 50kms^^, mean that the mea- 
surements should only provide a rough guide for the fainter 
galaxies. Further measurements can be found in the com- 
pilation of Bernardi et al. (2002), although the variety of 
sources from which these data were drawn mean that their 
reliability is less assured than in the other measurements. 

5.1 Results for individual galaxies 

In the following paragraphs we give a brief description of 
the observed kinematics for each individual galaxy in the 
sample. We discuss these result and compare them with the 
previous studies described above. 

5.1.1 NGC 1316 (FORNAX A) 

For this well-known merger remnant, Figure |21reveals gener- 
ally good agreement with the published results, although the 
new data reach to very much larger radii. The new data indi- 
cate a rather higher dispersion than the points at the largest 
radii of the previous measurements, but since these points 
were derived from the lowest S/N data in the older observa- 
tions, the conflict is probably not significant. The more ex- 
tended velocity dispersion profile from the new data for the 
first time reveals a seemingly-distinct feature of enhanced 
dispersion within the central 15 arcsec. In this galaxy, as 
might be expected in a merger, random motions dominate 
out to large radii in the new data, bringing into question 
whether this galaxy should really be classified as an SO sys- 
tem. 

5.1.2 NGC1380 

Figure |5| shows that the new data agree extremely well with 
the previous observations, but with much improved error 
bounds. The new data with their smaller errors and unfolded 
presentation give some indication of localized substructure 
that differs between the two sides of the galaxy. An enhance- 
ment in velocity dispersion within the central ~ 7 arcsec, 
first suggested by D'Onofrio et al. (1995), is strongly con- 
firmed with the new data. It is also notable from Figure IX2l 
that this region also shows structure in the skewness mea- 
sure, indicative of a distinct kinematic component. Interest- 
ingly, this feature coincides with photometric evidence for 
a distinct component which comprises a strong variation in 
isophotal ellipticity and position angle (Caon et al. 1994). 

5.1.3 NGC1381 

Figure |5| reveals a distinctive feature in this galaxy's rota- 
tional motion which was not apparent in the previous data. 
There is a correspondingly very strong signal in the /13 profile 
shown in Figure lX3l Such kinematic features are characteris- 
tic of a bar viewed edge-on (Bureau & Athanassoula 2005), 
so we would seem to have identified this edge-on galaxy, 
previously designated as unbarred, as a barred system. The 
new velocity dispersion profile diverges at large radii from 



the D'Onofrio et al. (1995) data, but this would seem to be 
because the latter has been artificially enhanced in regions 
where the dispersion approaches the resolution limit of the 
older data. 

5.1.4 IC1963 

Figure 121 shows very good agreement between new measure- 
ments of mean streaming velocity and those in D'Onofrio 
et al. (1995) for this rather nondescript SO system. The ve- 
locity dispersion measurements do not agree at all well, but 
once again this conflict can be attributed to the inadequate 
spectral resolution of the older data. 

5.1.5 NGC1375 

Once again. Figure |3 shows that the measured rotational 
motion for this galaxy agrees well with previous data, while 
the disagreement with earlier velocity dispersion data seems 
to arise from the limited spectral resolution of the older data. 
Support for this interpretation of the inconsistency comes 
from Kuntschner's (2000) value for the central dispersion of 
56 ± lOkms"^, in agreement with the new data. A further 
issue with the earlier analysis is revealed by Figure IA5I as 
the structure in and hi at small radii would indicate that 
there is more information in these dynamics than would be 
revealed by the earlier Gaussian fit. Once again, this struc- 
ture indicates the presence of a distinct kinematic compo- 
nent such as a bar. Photometric support for this conclusion 
comes from the work by Lorenz et al. (1993), who found 
that NGC 1375 possesses a complex structure characterised 
by boxy isophotes within a radius of 15 arcsec; such boxy 
isophotes are often associated with the kinematic signature 
of an edge-on bar (Kuijken & Merrifield 1995). 

5.1.6 NGC1380A 

The well-constrained measurements of rotational motion for 
this galaxy are fairly featureless, and, as Figurej^shows, do 
not reproduce the rather strange profile found by D'Onofrio 
et al. (1995). The similar absence of peculiarities in the 
and hi profiles shown in Figure lXSl suggest that the velocity 
structure in this galaxy is quite simple, so the smooth mean 
velocity curve is not unexpected. Once again, the new mea- 
surements of velocity dispersion lie systematically well below 
those in D'Onofrio et al. (1995), presumably due to the spec- 
tral resolution of the latter; Kuntschner's (2000) published 
central dispersion of 55±9 kms~^ fits rather better with the 
new measurements, but may also be running into resolution 
issues. 

5.1.7 ESO359-G002 

As Figure Oshows, this fainter lenticular galaxy has a slowly 
rising mean velocity curve, in good agreement with the pre- 
vious data from Graham et al. (1998) and despite being be- 
low our estimated spectral resolution. The previous velocity 
dispersion measurements are in reasonable agreement with 
the new data, but seem to lie systematically lower. In par- 
ticular, the sharp central peak in velocity dispersion, which 
is accompanies by a similar structure in hi (see Figure IXtIi . 
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Figure 2. Comparison between tlie new results for mean velocity, VloSi a^nd velocity dispersion, cr, and kinematics from the literature 
for NGC 1316, NGC 1380, NGC 1381 and IC 1963. The new results are shown as a continuous line with Itr errors as short dashed lines. 
The D'Onofrio et al. (1995) data are shown as triangles, the Longhetti et al. (1998) data are shown as stars, and the Graham et al. 
(1998) data are shown as squares. Horizontal dashed lines in the velocity dispersion plots indicate the resolution limit of the archival 
data. 



was not seen in the earlier data. In this case both, the veloc- 
ity curve and the dispersion profile should be viewed with 
some caution. 



58 ± 8kms"^ from Kuntschner (2000) and 46 ± 3kms"^ 
(Bernadi et al. 2002) are in accord with the new dispersion 
measurements. 



5.1.8 ESO358-G006 

There are no previous studies of the extended kinematics 
for this faint galaxy. Figure IA8I reveals generic featureless 
rotational motions that rise to a plateau at large radii, and 
a rather strange velocity dispersion profile with a slight cen- 
tral depression and a hint of a rising profile at large radii. 
The existing values for the central velocity dispersion of 



5.1.9 ESO358-G059 

The extended kinematics of this faint galaxy have not previ- 
ously been studied. Figure lXol shows the rather strange kine- 
matics for this system. It has a relatively well-behaved mean 
velocity profile that rises to an amplitude of ~ 50kms~^ at 
~ 5 arcsec, then declines slightly. The dispersion profile is re- 
markably fiat in the region where the rotation increases, then 
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Figure 3. As for Figure^ for NGC 1375, NGC 1380A and ESO359-G002 



starts to rise when the rotational velocities decline. The pub- 
lished values for the central velocity dispersion, 54 ±9kms^^ 
(Kuntschner 2000) and 46±3kms-^ (Bernadi et al. 2002), 
are in reasonable agreement with the new data. One possi- 
ble explanation for the enhanced dispersion and depressed 
rotation at larger radii would be the presence of a counter- 
rotating stellar disk (Rubin et al. 1992), but for such a faint 
galaxy we clearly do not have data of sufficient quality to 
test this hypothesis rigorously. 



6 CIRCULAR VELOCITY CALCULATION 

In order to go beyond the kinematic measurements of Sec- 
tion|K|to place these observations in their astrophysical con- 
text, we have to translate them into the related intrinsic 



dynamical properties of the galaxies. In particular, since 
we ultimately want to explore the TuUy-Fisher relation for 
these galaxies, we need to derive the circular velocity as a 
function of radius, Vc{R)- Although this quantity is clearly 
connected to the observed mean velocity profile, Vhos{R), 
the transformation from one to the other is not trivial, so 
we describe the procedure we have adopted in some detail. 

The first ingredient we need is the inclination, i, of the 
galaxy to the line of sight, in order to determine the fraction 
of rotational motion that will have been projected into the 
observed line-of-sight mean velocity. In principle, the incli- 
nation can be estimated by simply measuring the flattening 
of isophotes at large radii, but this approach is not advis- 
able as it uses only the lowest S/N photometry. Instead, we 
have carried out a full two-dimensional fit to images of the 
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galaxies using the GIM2D software (Simard et al. 2002) . We 
have applied this software to the publicly- available 2MASS 
Jf-band images of the objects in our sample (Jarrett et al. 
2003). The uniform quality of these data means that we 
can study all of the galaxies in a consistent manner, and 
the use of infrared light means that we are tracing the bulk 
of the population in these systems, avoiding any problems 
that might arise from localized star formation or dust ob- 
scuration. In addition to obtaining a best estimate for the 
galaxies' inclinations, this analysis also returns estimates for 
the parameters of a disk-plus-bulge model, including a bulge 
effective radius. Re, and a disk scalelength, -Rcxp. As we will 
see below, these photometric parameters are also useful in 
determining the properties of the rotation curve, so their 
values are presented in Table |H] 

To convert line-of-sight velocities into rotational mo- 
tions, we follow the approach developed by Neistein et al. 
(1999, hereafter N99). First of all, allowance has to be made 
for the extra degree of "roundness" contributed by the fact 
that the disks of SOs are not infinitely thin; as in N99, we 
assume that all galaxies would have an edge-on axis ratio of 
0.22, although this assumption turns out not to be critical 
to the results. 

The second major effect that must be corrected for is 
that arising from projection along the line-of-sight. Although 
stars close to the minimum radius along a line through a 
highly-inclined disk will have most of their rotational veloc- 
ity projected along the line of sight, those further from the 
centre of the galaxy will have most of their motion trans- 
verse to the line of sight, so their measured component of 
velocity will be small. This tail of low-velocity stars means 
that the measured line-of-sight mean velocity will lie sys- 
tematically below the rotational velocity of the stars. We 
have therefore applied a correction for this effect, which is 
described in Merrett et al. (2006). This technique can only 
be reliably applied in regions where the disk dominates the 
light, so we have used the photometric parameters derived 
above to ascertain the minimum radius, -Rlim, at which this 
correction can be made; the values for -Rlim are given in Ta- 
ble |3 Clearly, if -Rlim can't be estimated within the radius 
to which we have obtained data, then we will not be able 
to apply this correction, so will be unable to obtain a reli- 
able rotation curve. Where it can be applied, this correction 
amounts to ~ 20%; although it is clearly significant, it is 
still sufficiently small that its exact value is not critical. 

The last correction that we must apply is the difference 
between mean streaming velocity, V^, and circular rotation 
speed that arises because the stars have a random compo- 
nent to their motions as well as a rotational motion. This 
"asymmetric drift" can be corrected using the appropriate 
Jeans equation (Binney & Tremaine 1987), and here we fol- 
low the N99 implementation of this correction, by which 



2 



R 



Rc: 



- 1 



(2) 



where (jflt is the velocity dispersion at radius -R determined 
from a fit to the velocity dispersion profile using a low-order 
polynomial (see Appendix IXl. 

Strictly speaking, this low-order correction is only valid 
for rather cold disks with S> (Jfit- Some experimenting 
with simulated data revealed that the corrections remained 
fairly reliable for any galaxy with V^t/aat > 2, but that ar- 



tifacts from the correction would begin to appear in the ro- 
tation curve below this limit. We can meet this requirement 
for five of the sample galaxies, but for two (ESO 356-G006 
and ESO 358-G059), we must relax the limit to V^/aat > 1. 
In these cases, the asymmetric drift correction of equation|5| 
tends to be too large, so the resulting rotation velocities 
should strictly be taken as only an upper limit. In two cases 
(NGC 1316 and ESO 359-G002), random motions dominate 
at all radii, so this correction cannot even be attempted. In 
consequence, these two galaxies have to be excluded from 
subsequent analysis. 

After carrying out the above procedures, we can con- 
vert the measured mean line-of-sight velocity data points 
into discrete estimates of the local circular speed at the cor- 
responding radii. These estimates have two limitations: first, 
each point will have a significant error bar from all the prop- 
agated uncertainties; and second, the data will lie over only 
a limited range in radii since the rotation velocities cannot 
be determined in this way at small radii where random mo- 
tions dominate, and the galaxies become too faint at large 
radii to have their kinematics measured. We therefore can- 
not simply read off dynamical parameters such as the maxi- 
mum velocity, Knax, from these data. Instead, we effectively 
interpolate and extrapolate the limited data by modelling 
the contributions to the rotation curve by disk and bulge. 
Using the photometric parameters derived above, we model 
the disk component of the rotation curve using the formula 
for a pure exponential. 



v;f{R) = 4^GEoi?e.py1/o(j/)7fo(y) - hiy)K,iy)], 



(3) 



where G is the gravitational constant. So is the central sur- 
face density, 7o, Kq, Ii and Ki are modified Bessel functions, 
and y — -R/-Rcxp (Freeman 1970). For the bulge, we use the 
rotation curve of Hernquist (1990), 



Ige 



R^ 



R [R + aY 



(4) 



where Afbuigc is the mass of this component, and the scale- 
length a can be approximately related to -Rc by a/R^ — 1.82. 
The normalizations of these two components (effectively, the 
mass-to-light ratios of the photometric components) were 
then varied such that their sum produced the best fit to the 
circular speed data over the range that it is available. The 
result of this process is presented in Figure 0| we have not 
attempted this procedure for the two galaxies in which the 
circular velocities are very uncertain because the asymmet- 
ric drift corrections were too large. 

We can now use these fits to read off quantities of in- 
terest. For example, for a TuUy-Fisher analysis, we need 
the maximum rotation velocity Knax- The resulting values 
are presented in Tabled] for the two galax;ies for which we 
were unable to carry out the complete fit, we simply quote 
the maximum circular speed as derived from the data. Note 
that although we have generally derived Vmax from the fit- 
ting procedure, it is apparent from Figure |4]that the values 
returned do not differ greatly from the measured rotation 
speeds. This similarity arises because the new VLT obser- 
vations go deep enough to enable us, for the first time, to 
measure velocities in the flat part of the rotation curve. By 
reaching these radii, we enter a regime in which the maxi- 
mum velocity is reasonably tightly constrained irrespective 
of the details of the model fit. 
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Table 3. Structural and other important parameters of SO galaxies in the Fornax Cluster. 



Name 


i 




B/T 


Re 


.Rexp 


Sersic 




-ftout 




Mb 




["] 








1 


['■■ 


1 


n 


["] 


["] {[Re]} 


[kms~^] 








(2) 


(3) 


(4) 


(5) 


(6) 


(') 


(°) 




(10) 


NGC 1316 


43.2 


43.8 
42.4 


0.58 


0.62 
0.47 


36.0 


41.7 
24.8 


49.4 


52.7 
47.2 


2.9 


3.2 
2.3 




148.5 {4.1} 




-22.3 


(0.1) 


NGC 1380 


66.9 


67.6 
66.6 


0.58 


0.59 
0.57 


17.6 


18.0 
17.3 


36.3 


37.3 
35.9 


3.3 


3.4 
3.3 


35.7 


88.6 {5.0} 


309.6 (26.6) 


-20.6 


(0.1) 




82.5 


82.9 
82.2 


0.57 


0.58 
0.55 


7.4 


7.7 
7.3 


20.8 


21.6 
20.3 


3.1 


3.2 
3.0 


OO.Z 


ftA O /Q '71 
D4.Z \^C>. 1 } 




-19.0 


(0.1) 


IC 1963 


84.6 


84.7 
84.3 


0.39 


0.42 
0.36 


22.6 


24.9 
21.2 


15.8 


16.1 
15.5 


3.6 


3.7 
3.5 


5.5 


47.8 {2.1} 


164.7 (16.8) 


-18.5 


(0.1) 


NGC 1375 


67.5 


68.3 
66.3 


0.19 


0.21 
0.17 


3.2 


3.4 
2.9 


14.8 


15.2 
14.0 


2.2 


2.3 
2.1 


10.7 


35.9 {11.2} 


113.4 (19.6) 


-18.2 


(0.2) 


NGC 1380A 


78.1 


78.8 
77.4 


0.18 


0.24 
0.14 


10.9 


14.5 
8.1 


19.8 


20.4 
19.2 


3.7 


3.8 
3.6 


5.5 


47.9 {4.4} 


119.7 (3.0) 


-18.1 


(0.2) 


ESQ 358-G006 


65.9 


67.2 
64.8 


0.02 


0.04 
0.01 


1.5 


2.9 
0.2 


9.2 


9.6 
8.8 


3.0 


3.3 
2.7 


0.00627 


29.4 {19.6} 


155.8: (26.0) 


-17.5 


(0.2) 


ESQ 358-G059 


61.9 


63.9 
59.8 


0.72 


0.76 
0.65 


16.2 


18.4 
14.7 


2.3 


2.5 
2.2 


2.0 


2.2 
1.6 




14.2 {0.9} 


108.5: (53.8) 


-17.4 


(0.2) 


ESQ 359-G002 


49.4 


51.7 
47.2 


0.070 


0.14 
0.029 


36.0 


41.7 
24.8 


7.8 


8.2 
7.4 


3.1 


3.4 
2.9 


0.4 


13.9 {0.4} 




-17.3 


(0.2) 



Note: For all pertinent calculations, Hq = 70kms^-'" Mpc^"*". From (2) to (6), 99% confidence intervals arc presented; from (9) to (10), 

la errors between "()". Col (2), inclination angle with respect to the line-of-sight; col (3), bulge to total fraction; col (4), effective 
radius; col (5), exponential disk scale length; col (6), Sersic index; col (7), minimum radius where the "projected velocities correction" 
is applied; col (8), distance of the outermost bin in arc-seconds and in units of the effective radius; col (9), maximum circular velocity, 
values followed by ':' should be taken as an upper limit only; col (10), absolute magnitude in B-band, assuming a redshift of 0.0043 for 

the cluster, according to Madore et al. (1999). 



7 SUMMARY 



In this paper, we have presented kinematics derived from 
VLT long-slit spectral observations of 9 SO galaxies in 
the Fornax Cluster. The derived kinematic parameters axe 
the mean line-of-sight velocity, velocity dispersion, and the 
higher- moment /13 and /14 coefficients. The deep nature of 
the spectra mean that we have obtained these parameters 
out to distances equal or larger than 2 Re of the bulge for 
the majority of our objects, so we arc typically probing the 
disk-dominated part of the kinematics, allowing us to study 
both disk and bulge dynamics in some detail. 



A first qualitative pass through these kinematic data 
for the individual galaxies indicates that some of the existing 
data, due to both their lower signal-to-noise ratio and poorer 
spectral resolution, are unreliable. It also reveals that there 
is quite a lot of complexity in these systems, suggesting that 
their evolutionary histories may be relatively complex. With 
the quality of these new data, we should be able to obtain 
significant new archaeological clues as to how they evolved 
to their current states. As a first step in this direction, wc 
have estimated one of the principal defining dynamical char- 
acteristics of these systems, their maximum rotation speeds. 
In the next paper, we will use this quantity to explore the 
evolutionary clues in these systems' Tully-Fisher relation, 
while in a subsequent paper wo will exploit the high quality 
of the spectra to investigate the further clues that may be 
provided by the composition of their stellar populations out 
to large radii. 
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APPENDIX A: LINE-OF-SIGHT KINEMATICS: 

This appendix shows plots of the major-axis kinematics for 
the sample of 9 Fornax SO galaxies. The values of mean 
velocity, Vlos, velocity dispersion, cr, and higher-moment 
coefficients, and hi are plotted as a function of radius, 
and the location of each galaxy's effective radius, Rc, is also 
marked as a dashed line. The values of Rc and the position 
angle of the slit are annotated on each plot. The solid lines 
show the fit to the dispersion profile, crfit(-R), adopted in 
Section |S| 



so galaxies in Fornax: data and kinematics 



NGC 1316 (Fornax A) 



P.A.= 50.2" (HE-SW) 
R. = 35.0 " 
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Figure Al. The major-axis kinematics of NGC 1316. 



NGC 1380 



P.A.= 3.6° (N-S) 
K. = 17.6 " 
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Figure A2. The major-axis kinematics of NGC 1380. 
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NGC 1381 
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Figure A3. The major-axis kinematics of NGC 1381. 



IC 1963 



P.A.= 82.6° (E-W) 
K. = 23.6 " 
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Figure A4. The major-axis kinematics of IC 1963. 



so galaxies in Fornax: data and kinematics 



NGC 1375 



P.A.= 89.7° (E-W) 
K. = 3.3 " 




-10 

Hadius [■■] 



10 £0 30 40 



Figure A5. The major-axis kinematics of NGC 1375. 



NGC isaoA 



P.A.= -3.50' (N-S) 
R. = 10. S " 
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Figure A6. The major-axis kinematics of NGC 1380A. 



A.G. Bedregal, A. Arag on- Salamanca, M.R. Merrifield and B. Milvang- Jensen 



ESO 359-G002 



P.A - 47.3' (NE-SW) 
R. = 36.0 " 
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Figure A7. The major-axis Icinematics of ESO 359-G002. 



ESO 358-G006 



P.A.= 32.0° (NE-SW) 
R. = 1.5 " 
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Figure A8. The major-axis kinematics of ESO 358-G006. 
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Figure A9. The major-axis Icinematics of ESO 358-G059. 



